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1. Introduction 

Neutrinos have the weakest interactions among all known elementary parti- 
cles. They have also the smallest mass among all known massive particles. 
These two properties, on one hand, make it difficult to study neutrino 
properties directly, in particular, to measure their mass in laboratories. On 
the other hand, the same properties make neutrinos very important cos- 
mologically and, at the present time, measuring neutrino masses looking 
at the sky seems more promising than terrestrial experiments. Significant 
cosmological role played by neutrinos arises from their large number den- 
sity. Neutrinos are the second most abundant particles in the universe, 
after photons in Cosmic Microwave Background Radiation (CMBR) with 
the number density n 7 w 410/cm 3 . According to the standard cosmol- 
ogy, the universe is filled, in addition to CMBR, by the Cosmic Neutrino 
Background Radiation (Cz^BR) with the present-day number density: 



56/cm 3 (1) 



for any neutrino flavor a = e,/j,,T. It is usually (but not always) assumed 
that neutrinos are not degenerate (i.e. their chemical potentials are zero or 
negligibly small) and the number densities of neutrinos and antineutrinos 
are equal. 



1 
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However, contrary to well observed CMBR, the existence of CjaBR is 
only a theoretical prediction, and though practically nobody has any doubts 
of that, direct observation of CVBR is still missing and it seems that there 
is no chance for direct registration of cosmic neutrinos today and maybe 
even in the foreseeable future. Thus, one has to rely on indirect methods 
studying features imprinted by cosmic neutrinos on: 

(1) formation and evolution of astronomical Large Scale Structure 
(LSS); 

(2) angular fluctuations of CMBR; 

(3) light element abundances created at Big Bang Nucleosynthesis 
(BBN); 

(4) propagation of Ultra High Energy Cosmic Rays (UHECR) . 

In what follows we will consider the first three subjects only. One can 
find a more detailed discussion of the issues presented below, as well as of 
the problem of interaction of UHECR with C^BR, and hopefully a complete 
list of references in the review 1 . 

2. Early History of Neutrinos. 

When the temperature of the cosmic plasma was above a few McV, neutri- 
nos were in thermal equilibrium with the electromagnetic component of the 
plasma, i.e. with photons and e + e _ -pairs. Electronic neutrinos decoupled 
from e + e~-pairs when the temperature dropped below w 1.9 MeV, 
while v ^ and v T decoupled a little earlier at w 3.1 McV. At the mo- 
ment of decoupling and later down toT« m e = 0.511 MeV temperatures 
of neutrino and electromagnetic components were equal, T v = T 7 . Below 
T = m e the annihilation of e+e~-pairs heats up photons and themselves, 
while leaves neutrino temperature intact. As a result the initial equilibrium 
ratio of neutrino-to-photon number densities becomes diluted by the factor 
4/11: 

n„ + n p _ 3 _3_ 

n 7 ~ 4 ^ 11 U 

From this result the present-day neutrino number density (1) is obtained. 
Any additional energy release after which might increase the photon 
number density would correspondingly diminish 

Neutrino spectrum is close to the equilibrium form: 

= [exp(p/T v -0 + If 1 (3) 
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where p is the neutrino momentum, T„ is the temperature, and £ = p,/T 
is dimcnsionless chemical potential; for adiabatic expansion £ remains con- 
stant. Usually chemical potentials of different neutrino species are assumed 
to be negligible, at the level of baryonic asymmetry, 1CP 9 , though much 
larger values, even close to 1, are not excluded. Moreover, there exist 
theoretical models which predict a small baryon asymmetry and simul- 
taneously large lepton ones. Equilibrium with respect to ^-annihilation 
enforces + £p = but if neutrino charge asymmetry is generated at low 
temperatures this condition may be violated. 

One more comment is in order: if neutrino mass is non-negligible in 
comparison with the temperature, the spectrum (3) is non-equilibrium be- 
cause the latter contains exp(E/T) but not exp(p/T). At the present time 
T v w 0.714T 7 w 1.68 • 10~ 4 cV. And if the neutrino mass is larger than this 
value the deviations from the usual Fermi distribution may be significant. 
This must be true at least for two out of three neutrinos because from the 
atmospheric neutrino anomaly Sm 2 = (2 — 5) • 10~ 3 eV 2 and from the solar 
neutrino deficit and KAMLAND data Sm 2 = (6 - 7) • 1CT 5 cV 2 . 

However, even at high temperatures (T 3> m v ) a deviation from equi- 
librium was non-negligible. Indeed, the decoupling of neutrinos from e+e~- 
pairs is not instantaneous and the annihilation e+e~ — > vv at T < m e 
would heat the neutrino component of the plasma and distort its spectrum. 
According to analytical estimates of rcf . 2 the spectral distortion for v e has 
the form: 



Most accurate numerical solution of the kinetic equation that governs non- 
equilibrium corrections to the neutrino spectrum was performed in rcf. 3 . 
According to the calculations the excess of energy density of v e and v^ T 
are respectively 5p Ve / ' p v = 0.9% and 5p v Vr / p v = 0.4%. Together with the 
plasma corrections which diminish the energy density of the electromag- 
netic component, the total relative rise of neutrino energy density reaches 
approximately 4%. This phenomenon has very little impact on production 
of primordial 4 Hc but can be observable in the shape of the angular fluc- 
tuations of CMBR in the forthcoming Planck mission. If observed, then 
together with BBN, it would present evidence of physical processes which 
took place in the universe when she was about 1 sec old. The corresponding 
red-shift is about 10 10 . More details and references relevant to the subject 
of this section can be found in the review 1 . 
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3. Gerstein-Zeldovich and LSS Bounds on Neutrino Mass 

Knowing the present-day number density of relic neutrinos one can easily 
calculate their energy density and obtain an upper limit on their mass. 
Such bound was derived in 1966 by Gerstein and Zeldovich 5 The result 
was re-derived 6 years later by Cowsic and Mc Lelland 6 but in their work 
the effect was overestimated by the factor 22/3. In contemporary form the 
limits reads: 

J2 m ^ <94eV£Wi 2 (5) 

a 

where the sum is taken over all neutrino species, a = e,fi,r; ft = p/p c is 
the cosmological mass fraction of matter, p c — 10.5 h 2 keV/cm 3 , and h is 
the dimensionless Hubble parameter, h = H/100 km/sec/Mpc w 0.7. 

According to the different and independent pieces of astronomical data 
< 0.3 and correspondingly ^ a m„ a < 14 eV. Since the data on neutrino 
oscillations 7 show that the neutrino mass difference is much smaller than 
eV, the mass of any neutrino flavor should be below 4.7 eV. 

This limit can be further improved if one takes into account a possible 
role that massive neutrinos might play in the process of formation of the 
large scale structure of the universe. The point is that in neutrino dom- 
inated universe all the structures on the scales smaller than the neutrino 
free streaming length, l v , should be erased. The latter can be estimated as 
the distance that neutrinos could free-stream in the universe before they 
became non-rclativistic. The mass inside the free-streaming volume can 
be estimated as M v w 5 • 10 17 M Q (leV/m^) 2 . If density perturbations are 
common for all particles (this is true for adiabatic perturbations created 
by inflation), then neutrino out-stream would leave behind less power at 
small scales inhibiting structure formation at these scales. The larger is Q u 
the larger is the effect. Moreover, the larger is the fraction of relativistic 
(hot) dark matter the later structure formation begins. Hence observa- 
tion of structures at large red-shifts allows to conclude that f2„ < 0.1 and 
m v < 1.5 eV 8 . The analysis of the recent data from 2dF Galaxy Redshift 
Survey 10 permitted to put the limit m„ a < 2.2 eV or individual masses 
should be below 0.73 eV. As argued in ref. 9 detailed analysis of structure 
formation by the Sloan Digital Sky Survey would be sensitive to neutrino 
mass at the level of (a few) x 0.1 eV. 

This limits are based on certain assumptions about the form of the 
spectrum of perturbations, their character (adiabatic or isocurvature), and, 
what is probably the safest one, about neutrino interactions. A question 
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may arise in this connection whether neutrinos can supply all the dark 
matter in the universe if we relax any or all these assumptions, introducing 
an arbitrary shape of the perturbation spectrum and/or new neutrino in- 
teractions (of course inside the established limits). The answer would still 
be negative because of the Tremain-Gunn limit 11 . This limit manifests 
quantum mechanical Fermi exclusion principle at the kiloparsec scales and 
demands that neutrinos must be heavier than roughly speaking 100 eV if 
they form all dark matter in galaxies. So it seems that the only way out 
is to make a crazy assumption that neutrinos are bosons and not fermions 
which in light of the discussed today search of CPT violation maybe looks 
not so crazy because CPT-theorem is heavily based on the standard relation 
between spin and statistics. 

4. Neutrino Mass and CMBR 

The impact of neutrinos on the shape of angular spectrum of CMBR tem- 
perature is based on the following two effects (for details see e.g. the recent 
review 12 ). First, a change of the energy density of relativistic matter 
(neutrinos) would change the cosmological expansion regime and this in 
turn would change the physical size of the horizon at recombination. It 
would shift the positions of the acoustic peaks in the temperature fluctua- 
tions. This effect is relatively weak. More important is another one which 
changes the heights of the peaks because an admixture of relativistic matter 
gives rise to an amplification of acoustic oscillations since relativistic matter 
creates gravitational force which varies with time creating resonance ampli- 
fication of the acoustic oscillations. These phenomena permit to measure 
neutrino mass and the number of neutrino families (or in other words, the 
energy density of neutrinos) in the recombination epoch. 

The recent measurements of the angular fluctuations of CMBR by 
WMAP 13 , together with the analysis of LSS by 2dF, permitted to im- 
pose surprisingly strong upper limit on neutrino mass 



that is for mass degenerate neutrinos m v < 0.23 eV at 95% CL; see also 
the discussion of this result and of the role of primers in ref. 14 . 

Immediately after publication of the first WMAP data, there appeared 
several papers 15 where the number of neutrino families were evaluated. 
Their results are, roughly speaking, N u = 1 — 7 depending upon the analysis 




(6) 



a 
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and priors. Anyhow N v > and thus CMBR confirms, independently on 
BBN, that Cz^BR (or some other relativistic background) indeed exists. 

The accuracy of the data is not good enough to compete with the deter- 
mination of N u from BBN (see sec. 5) to say nothing about a registration 
of 4% addition to the neutrino energy discussed in sec. 2. However the 
forthcoming Planck mission may be sensitive to this effect and overtake 
BBN in determination of N u . One should keep in mind that the impact of 
neutrino energy on BBN and CMBR may be different depending upon the 
form of the spectral distortion of electronic neutrinos f Ve {E). The effect, 
which was calculated for massless neutrinos, also depends upon the value 
of neutrino mass. This may be essential at the recombination epoch. 

5. Neutrinos and BBN 

Abundances of light elements ( 2 H, 3 He, 4 He, and 7 Li) produced at BBN at 
T = 1 — 0.07 MeV (t=l-200 sec) depend upon the following quantities: 

(1) Number density of baryons, 7710 = 10 10 ns/n 7 . In the previous 
century the value of this parameter was determined from BBN it- 
self through comparison of the predicted deuterium abundance with 
observations. Now after measurements of CMBR angular fluctua- 
tions by MAXIMA, BOOMERANG and DASI 16 and confirmed by 
WMAP 13 this parameter is independently fixed at 7710 = 6 ± 0.3. 

(2) Weak interaction rate which is expressed through the neutron life- 
time, now well established, 885.7 ± 0.8 sec 17 . 

(3) Cosmological energy density at the period of BBN. The latter is 
usually parametrized as the effective number of additional neutrino 
species AN U = N v — 3. This type of parametrization is flawless 
for relativistic energy, while for another form of energy (e.g. non- 
relativistic or vacuum-like) its effect on the production of different 
light elements may deviate from that induced by neutrinos. 

(4) Neutrino degeneracy, i.e. a possible non-vanishing values of neu- 
trino chemical potentials £ . While for and v T non-zero £^ )T are 
equivalent to an increase of N Vl degeneracy of electronic neutrinos 
has a much stronger impact on BBN because the frozen neutron-to- 
proton ratio is exponentially sensitive to the magnitude (and sign) 
of chemical potential of electronic neutrinos, n n /n p ~ exp(— £„ e ). 

(5) Energy spectrum of v e . If the spectrum is distorted then the equi- 
librium of reactions nv e <-» and ne + <-> pi7 e is shifted and could 
result either in an increase or decrease of n n /n p depending upon 
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the form of distortion. 

The upper bound on the effective number of neutrino families, found 
from BBN, had a rather strong time evolution and it changed from 
AiV„ = 3 — 5 in the original papers 18 down to 0.3 and even to 0.1 in 
subsequent literature. The existing data on the light clement abundances 
is still controversial 19 and possibly the situation described in rcf. 20 still 
remains true, which is roughly the following n v — 3 ± 0.5. 

The most accurate upper bound on chemical potential of a single neu- 
trino species under the conservative assumptions that A7V < 1 and other 
chemical potentials are vanishingly small was obtained in ref. 21 and reads: 
|£i/„,i/ T | < 1-5 and |£„J < 0.1. This bound would be relaxed if a conspir- 
acy between different chemical potentials is allowed such that the effect of 
a large \£ v ;I/ J is compensated by |£„J (or vice versa). The best limit is 
presented in ref. 22 where the combined inputs from BBN and CMBR have 
been used: I^^J < 2.6 and |£„J < 0.2. 

The bounds quoted above are valid if neutrinos are not mixed. Other- 
wise, charge asymmetry of a certain neutrino flavor would be redistributed 
between all neutrino species and the bound would be determined by the 
most sensitive asymmetry of v e . This problem was analyzed in ref. 23 
numerically and analytically in ref. 1 , where it was shown that for Large 
Mixing Angle solution of the solar neutrino deficit the transformation is 
quite efficient and equilibration of all chemical potentials by oscillations 
is achieved. This leads to the common bound for all neutrino chemical 
potentials 

16*1 < 0.07 (7) 

Similar investigation both analytical and numerical was also performed in 
the papers 24 . 

If a new sterile neutrino (or several sterile species) mixed with ac- 
tive ones exists, the impact of neutrino oscillation on BBN would be 
more interesting than in the case of active-active mixing. The v s <-> v a - 
transformations would excite additional neutrino species leading to an in- 
crease of the effective number of neutrinos, AN > 25 > 26 , could distort 
the spectrum of u e 26 ' 27 and could generate a large lepton asymmetry in 
the sector of active neutrinos by MSW-resonance 28 . (More references and 
discussion can be found in .) 

In the non-resonance case, i.e. for m„ s > m Va it is rather easy to esti- 
mate the production rate of sterile neutrinos in the early universe through 
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oscillations and to obtain the following bounds on the oscillation parame- 
ters 29 : 

(Sml^/eV 2 ) sin 4 20^ = 3.16 • 1CT 5 [log (1 - AN,)} 2 (8) 
(Sml^JeV 2 ) sin 4 2ff£? = 1.74 • 1(T 5 [log (1 - AA^)] 2 (9) 

In this result a possible deficit of v e created by the transformation of the 
latter into v s , when refilling of v e by e + e~-annihilation is already weak, is 
not taken into account. This effect would strengthen the bound. 

Numerical solution of kinetic equations governing neutrino oscillations 
in the early universe was performed in rcf. 30 under assumption of kinetic 
equilibrium, so the neutrinos are described by a single momentum state with 
the thermally average value of the energy E — 3.15T. For the non-resonant 
case the obtained results are somewhat stronger than those presented above, 
(8,9). For the resonant case it is questionable if thermal averaging is a 
good approximation because the position of resonance depends upon the 
neutrino momentum. One more complication is that now is known that all 
active neutrinos are strongly mixed and their mutual transformation should 
be taken into account together with [v a — v s ) -oscillations 31 . The active 
neutrino mixing noticeably changes the previously obtained cosmological 
limits for mixing with v s . Moreover, the solution of momentum dependent 
kinetic equations shows that kinetic equilibrium is strongly broken (at least 
for some values of the oscillation parameters) and, in particular, the spec- 
trum of v e is distorted leading to a shift of the n/p-ratio for interesting 
values of mass difference between v a and active neutrinos. The calculations 
are complicated by the appearance of one, two or even three resonances if 
sterile neutrino is lighter than one, two or all three active ones. Anyhow in 
the resonance case the cosmological bounds on the mixing between sterile 
and active neutrinos are considerably stronger than those in non-resonance 
case. Thus if a large mixing to v s is discovered it would mean that the 
lcpton asymmetry of the universe is non-negligible 31 , because the latter 
might "cure" the effect of v a — v s oscillations on BBN. 

6. Conclusion 

Thus we see that cosmology is becoming sensitive to the values of neutrino 
masses approaching V 6m 2 . So one may hope that neutrino will be the first 
particle whose mass will be measured by astronomers by the combined data 
from CMBR and LSS. 

The number of additional neutrinos at BBN is limited by 0.5 (though 1 
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is probably still not excluded) with a possibility to improve this limit down 
to 0.1. 

The observed in experiment strong mixing between active neutrinos 
allows cosmological lepton asymmetry to be relatively low, smaller than 

0. 1. This excludes, in particular, cosmological models where large chemical 
potential of neutrinos might be essential for large scale structure formation. 

A possible mixing between active and sterile neutrinos is restricted by 
BBN at much stronger level than by direct experiment. 

Planck mission may be sensitive to additional contribution to neutrino 
energy density at per cent level and thus will be able to trace physical 
processes in the universe at red shift of 10 10 . 

I am grateful to F. Villantc for critical comments. 
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